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Abstract 


A  study  of  the  effects  of  cell  temperature  and  pressure  on  anode  polarization  was 
achieved  in  this  paper.  A  Langmuir  and  Temkin  isothenn  was  used  to  study  the  CO  adsorption 
equilibrium  and  the  Heyrowski-Volmer  mechanism  for  fb  oxidation  on  a  Pt/Ru-C  catalyst.  The 
anode  overpotential  curves  were  separated  into  four  regions  so  that  the  adsorption  rate  constant 
for  H2  on  Pt/Ru-C  alloy  could  be  determined.  This  separation  allowed  for  the  generation  of  a  CO 
isotherm  and  it  provided  Lagmuir  parameters  to  describe  the  CO  adsorption  equilibrium  as  a 
function  of  temperature  and  pressure. 


Introduction 

Recently  Murthy  et  a  I.1  show  dramatic  changes  is  steady  state  performance  of  a  PEMFC 
when  the  temperature  was  changed  from  70°C  to  90°C  with  500  ppm  CO  in  TL  (CO/H2).  They 
also  show  substantial  changes  in  the  rates  of  poisoning  and  recovery  for  this  20°C  temperature 
change.  These  data  were  obtained  using  a  Pt/Ru-C  alloy  electrode.  Also  recently  Springer  et  al. 
showed  equations  to  study  both  steady  and  transient  coverage  of  CO  on  membrane  and  electrode 
assemblies  (MEA)  with  platinum  catalysts.  These  equations  depend  on  a  large  number  of 
thermodynamic  and  kinetic  parameters  for  adsorption  and  electrochemical  oxidation  of  both 
hydrogen  and  CO.  They  gave  a  set  of  parameters  that  fit  their  data  and  with  these  parameters 
and  their  equations  they  were  able  to  obtain  agreement  with  steady  state  voltage-current  data  as  a 
function  of  CO  concentration  and  system  pressure.  However  they  did  not  investigate  the  effect 
of  temperature  and  how  these  parameters  differ  for  an  alloy,  particularly  Pt/Ru-C.  Other  studies 
have  shown  that  there  is  a  significant  change  in  cell  performance  due  to  a  change  in  temperature 
at  low  CO  concentrations  (i.e.,  20  ppm)  .  There  have  also  been  several  studies  showing  the 
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effect  of  an  alloy,  such  as  Pt/Ru-C  on  the  overall  poisoning  effect  of  CO  in  the  PEMF  '  .  In  this 
paper  we  focus  on  the  change  in  these  parameters  as  a  function  of  temperature  and  pressure  for 
both  steady  state  and  transient  behavior  during  exposure  to  CO  in  PC. 

We  are  motivated  by  recent  anodic  polarization  data  that  were  reported  for  relatively  high 
concentrations  of  CO/H2  on  a  PRIME  A  MEA  Series  551’8’9.  Current  work  in  our  laboratory  is 
directed  at  the  effect  of  temperature  and  pressure  on  the  steady  state  and  the  rates  of  recovery 
and  poisoning  during  exposure  to  concentrations  of  500,  3000  and  10000  ppm  CO/H2. 
Specifically  data  were  reported  at  202  kPa  for  70°C  and  90°C  in  addition  to  the  steady  state 
results1.  Rates  of  poisoning  and  recovery  for  step  changes  in  the  CO  levels  in  hydrogen  were 
also  reported.  Examples  of  these  transient  data  are  shown  in  Figure  1  where  the  cell  was 
operated  at  a  fixed  current  density  of  600  mA/cm  with  an  initial  stream  of  pure  hydrogen.  The 
decrease  and  increase  in  voltage  resulted  from  changes  in  the  CO  concentration  in  the  anode 
feed.  We  also  have  polarization  curves  for  concentrations  of  500  and  3000  CO/FE  at  101  kPa 
and  70°C8.  However  since  major  emphasis  was  on  providing  poisoning  and  recovery  rates  rather 
than  for  studying  the  effect  of  cell  temperature  on  cell  performance,  the  framework  provide  for 
that  work  is  limited  and  will  be  given  a  lower  priority  in  estimating  the  parameters. 

The  objective  of  this  study  is  to  provide  parameters  for  the  equations  that  describe  a 
Pt/Ru-C  catalyst  used  in  high  performance  MEA  such  as  PRIMEA  MEA  Series  55.  These 
parameters  will  allow  the  description  of  performance  as  a  function  of  temperature  and  pressure. 
The  sensitivity  of  the  predictions  to  changes  in  the  parameters  will  also  be  presented. 
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Theory 


Giorgi  et  al6  recently  proposed  the  Heyrovski-Volmer  mechanism  for  the  oxidation  of  H2 
on  a  Pt/Ru-C  catalyst.  Thus  a  set  of  reactions  for  the  alloy  should  be  written: 


kfc 


C0+PfO(Pf-C0j 

k  fcbfc 

[1] 

H2  +  Pt<2,(Pt-H2J 

kfcbfl, 

[2] 

{Pt-HlaJ,)<P>2H'+2e~+Pt 

[3] 

kec 

H20  +  [Pt  -  COads  )<^>  C02  +  2 H+  +  2e~  +  Pt 

[4] 

Then  the  rate  expressions  for  the  coverage  of  CO  and  H2  on  the  Pt  sites  of  the  alloy  can 
be  written: 


Where: 
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vah  =  n, 


for  i 


total 


\_ 

2 


(Va+Uco) 


for  hotel  =  C  +  L  [12] 


This  mechanism  is  in  contrast  to  the  Tafel-Volmer  mechanism  of  fh  oxidation  on  Pt  used  by 
Springer  et  al.  .  Note  that  in  Reaction  1-4  and  Equations  5  and  6,  we  have  followed  the  notation 
of  Springer  et  al.  and  thus  one  can  recognize  that  the  first  and  second  terms  in  Equation  6  do  not 
include  a  power  ‘n’  which  should  be  equal  to  2  for  the  Tafel-Volmer  mechanism.  The  molar 
area  density  of  catalyst  sites,  ct,  has  units  of  mol/cm2.  The  adsorption  rate  constants  kfc  and  kjh 
are  for  CO  and  TE  adsorption  respectively  with  units  of  mol/s-cnr-atm.  bfc  ( l/Keqxo)  and  bfh 
( 1/Keqj, )  represent  the  back-to-forward  CO  and  Ho  adsorption  ratio  with  Pa  units.  Also  kec  and  keir 
are  the  electrochemical  oxidation  rate  constants  for  CO  and  H2  oxidation  respectively,  with  units 
of  mol/cm  -s. 

Equations  7  and  8  are  the  current  density  expressions  for  CO  and  H2  as  a  function  of  their 
respective  overpotentials  and  coverage  and  Equations  9  and  10  are  their  respective  exchange 
current  density  expressed  in  terms  for  there  respective  equilibrium  coverage.  Equations  1 1  and 
12  are  the  overpotentials  of  CO  and  H2  in  terms  of  the  experimental  anode  overpotential 
assuming  that  the  equilibrium  potential  for  H2  oxidation  is  zero  and  that  the  anode  overpotential 
for  H2  oxidation  in  the  absence  of  CO  is  negligible.  At  the  exchange  current  density,  therefore  in 
absence  of  H2  and  CO  oxidation  (i.e.,  at  i,otai  =  0)  the  following  expressions  can  be  derived  for 
the  steady  state  equilibrium  Oco  and  Oh'. 
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Note  that  these  are  their  equilibrium  coverage.  The  sum  of  the  individual  current  densities,  h 
and  ico,  is  the  total  or  system  current. 


CO  Isotherm  Models 

The  two  most  referred  to  models  for  CO  adsorption  onto  Pt  in  a  Pt-C  or  Pt  alloy  catalysis 
are  the  Langmuir  and  Temkin  isothenns.  If  we  take  Equation  13  in  the  absence  of  TL  we  get  the 
well  known  Langmuir  model,  shown  below. 


0, 


CO,  pure 


[14] 


'fc 


CO 


+  1 


13 

Taking  the  relationship  presented  by  Chafik  et  al.  between  the  equilibrium  adsorption 
coefficient  ( Keq=l/bfc )  and  the  adsorption  energy,  bfc  is  defined  as  the  following: 


,  1 
bfc  =^eXP 


RT 


[15] 


Where  A  is  the  pre-exponential  tenn  in  their  expression  and  Eads  is  the  adsorption  energy.  Note 
that  as  in  the  work  by  Chafik  et  al.  ,  A  was  detennined  at  a  fixed  temperature  (600  K).  A 
Langmuir  isotherm  model  assumes  that  Eads  is  constant  with  Oco-  Therefore  Equation  1 5  can  be 
used  as  it  is  for  a  single  value  of  Eads.  However  in  the  case  of  Temkin  isotherm  it  is  assumed  that 
Eads  changes  linearly  with  Oco  according  to  the  following  expression: 

E*=E,  114] 


Where  Eo  is  Eads  at  Oco  =  0  and  a  is  the  fitting  parameter  that  is  independent  of  temperature.  If 
we  take  AE  to  be  the  difference  between  Eads  at  Oco  =  0  and  at  Oco  =  1  then  from  Equation  14  we 
get  AE  =  oEq.  Upon  substitution  into  Equation  15  the  following  expression  for  bfc  as  a  function 
of  Oco  is  obtained: 
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bfc  =  bfc0  exp 
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RT 
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RT 


[15] 


[16] 


Note  that  this  expression  is  similar  to  that  given  by  Springer  et  al?  Where  AE  is  SAGco  in 
notation  used  by  Springer  et  al.  .  Also  note  that  bfco  depends  on  temperatures,  which  was  fixed 
in  the  work  by  Springer  et  al?. 


Rate  of  Poisoning  and  Recovery 

All  poisoning  and  recovery  experiments  in  Ref.  1  and  8  where  conducted  at  a  fixed 
current  and  their  results  expressed  as  cell  voltage  as  a  function  of  time.  Hence  taking  the 
derivative  of  the  sum  of  the  two  current  densities  (Equation  7  and  8)  and  equating  them  to  zero, 
the  following  expression  can  be  derived  for  change  in  the  anode  overpotential: 


Where: 
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Equation  17  gives  an  expression  that  supports  the  assumption  that  the  poisoning  and  recovery 
rates  are  governed  by  surface  coverage  phenomena1, 8. 

The  anode  potential  curves  shown  in  Ref.  8  were  detennined  by  subtracting  the  cell 
voltage  in  the  presence  of  CO  from  that  of  neat  H2  at  a  fixed  current.  This  was  done  with  the 
assumption  that  the  mass  transfer  controlled  ovepotential  from  02  diffusion  is  only  a  function  of 
current  and  that  the  anode  overpotential  from  H2  oxidation  during  neat  H2  is  negligible.  The 
expression  below  describes  the  anode  overpotential  following  the  above  assumptions: 

V <,{*)=  Kell \h 2(CE)  “  Vcell\cO/H1  W  [24] 


Where  V, 


eell\H2(CE) 


is  the  cell  voltage  for  neat  hydrogen  at  a  particular  current  density  on  clean 


electrode  (CE).  Note  that  its  value  does  not  change  because  they  depend  on  current,  which  is 
fixed  during  the  experiment.  Therefore  the  change  is  brought  about  by  the  anode  polarization, 
which  depends  on  the  rate  of  coverage  of  H2  and  CO. 

After  some  poisoning  experiments,  Murthy  et  <://. 1,8  observed  that  the  cell  did  not  recover 
to  its  original  value  and  thus  we  write  Vcell\  ^  =  Vcell\  E) .  Where  Vcell\H  (pE)  corresponds  to 

the  cell  voltage  for  neat  hydrogen  on  a  poisoned  electrode  (PE).  We  speculate  that  some  sites 
were  non-recoverable  from  exposure  to  CO  and  thus  additional  polarization  can  be  written: 
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[25] 


77  =  V  -V 

la\  H2(PE)  cell\ H2(CE)  cell\ H2(PE) 

That  is,  Va\H  (PE)  's  the  anode  overpotential  resulting  from  the  FF  oxidation  reaction  for  an 
exposed  electrode.  This  value  changes  depending  on  the  degree  of  CO  poisoning. 


Results  and  Discussions 

These  equations  depend  on  the  twelve  parameters;  bfc  (8(AGco),  AGC0|  ,A)  A;,,  kfc, 

kfih  kec,  keh,  Pi,.  pc0,  Uco9  and  ct.  Eleven  of  parameters  were  estimated  using  the  steady  state  data 
from  and  the  remaining  from  transient  data  from  Ref.  1.  To  estimate  the  values  of  these  eleven 
parameters  each  steady  state  anode  overpotential  perfonnance  curves  was  separated  into  four 
distinct  regions  shown  in  Figure  2.  The  data  in  Figure  2  is  for  3000  ppm  CO  in  H2  at  70°C  and 
202  kPa  from  Ref.  1 .  The  sequence  for  these  estimation  perfonned  on  the  data  from  Ref.  1  for 
the  conditions  listed  in  Table  I  is  as  follows:  bn,  and  keh  are  obtained  from  Region  I;  S(AGco), 
AGCO|0  ()  and  kfh  are  obtained  next  from  Region  II;  finally  kfC,  kec  and  Pco  are  obtained  from 

Region  III.  For  three  of  the  four  remaining  points  A,  Uco9  and  Ph  their  values  are  assumed  in  a 
manner  consistent  with  the  data  of  Ref  1.  Finally  ct  is  obtained  from  the  transient  recovery  and 
poisoning. 

The  objective  function  show  below  was  used  in  estimating  the  parameters  described  in 
detail  below: 


Nc  Np  Nt  Nd 

<^=ZZZZ(' 

C-\  P= 1  T= 1  D= 1 


obs  pred  p 

J1  a  ~71  ‘ 


. 


D,T,C 


[26] 


Where  P,  C,  T  and  D  are  the  CO  concentrations,  cell  Pressure,  cell  temperature  and  data  range 
respectively. 
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Estimating  bjh  and  ke> Region  I  (Low  Polarization  Region) 

The  anode  overpotential  in  Region  I  shown  in  Figure  2  is  quite  low.  Which  indicates  that 
the  H2  oxidation  reaction  has  enough  available  free  surface  for  the  reaction  even  after  many  sites 
have  been  covered  with  CO2, 5.  In  this  region  we  assume  that  the  coverage  of  CO  remains  at 
equilibrium  throughout  Region  I  and  that  i  is  linear  with  r|a  at  low  current  density  values. 
Inherent  in  that  assumption  is  fact  that  0H  has  to  be  a  low  enough  for  0co  to  remain  at  its 
equilibrium  value  even  after  a  potential  is  applied.  Also  this  will  indicates  that  0h  is  also  low  in 
the  absence  of  any  competing  adsorbent.  Therefore  a  value  of  2%  was  assumed  for  0H*  which  is 
the  equilibrium  H2  coverage  in  the  absence  of  CO.  Then  considering  only  the  steady  adsorption 
of  FL  without  the  presence  of  CO,  from  Equation  11,  a  relationship  of  the  following  can  be 
formed  for  bfj,: 

[27] 

p„,  e» 

The  ratio  bjj/PH2  =10  for  the  assumed  value  of  0h*.  We  detennined  bn,  based  on  that  PH2  at  202 
kPa  since  for  which  more  data  was  available.  Note  that  bf 7,  was  considered  to  be  constant  with 
temperature  since  there  was  no  significant  difference  in  cell  perfonnance  under  neat  H2  operation 
at  70  and  90°C. 

Next  we  determined  keh  using  the  value  of  bn,  determined  from  above.  The  data 
corresponding  to  the  condition  of  500  ppm  CO/FL  at  90°C  shown  in  Table  II  was  used  in 
detennining  keh.  Assuming  that  ih  is  linear  with  T|a  and  that  0H  =  0H,eq  throughout  Region  I,  the 
linear  form  of  Equation  20  shown  below  was  used. 
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The  value  of  io,h  was  varied  such  that  a  minimum  OF  was  obtained  and  from  which  keh  was  then 
detennined  using  Equation  10. 


Estimating  b/c:  Region  II  (steady  state  total  rate  limited  by  available  surface) 

In  this  region  the  sites  covered  by  the  adsorbed  CO  limit  the  H2  oxidation  reaction  hence 
the  sharp  increase  in  the  anode  overpotential.  At  this  point,  the  available  free  sites  for  the 
oxidation  reaction  are  used  and  the  H2  surface  coverage  approaches  zero.  We  then  use  imtai  equal 
to  a  surface  limited  current  density,  iieco ■  When  this  occurs,  the  anodic  polarization  continues  to 
increase  without  bound  until  the  adsorbed  CO  is  removed  from  the  surface  by  oxidation. 
Therefore  Oco  is  constant  up  from  i]a  =  0  until  an  overpotential  corresponding  to  iieco ■  Applying 
these  limitations  to  the  steady  adsorption  forms  of  Equations  5  and  6  including  electrooxidation 
of  H2  and  solving  those  systems  of  equations  for  Oco  and  Oh  yeilds: 


he  I  b , 

l'°«>  =  k  P  ’ 

IF  *  H\  bfc+Pco 


Where  iieco  is  the  total  rate  fixed  by  the  non-oxidized  CO  on  the  surface  of  the  electrode.  Notice 
that  the  above  expression  is  similar  to  that  given  by  Springer  et  al.~  however  the  exponent  ‘n’ 
does  not  appear  on  the  bracketed  expression  due  to  the  previously  mentioned  assumption  for  FE 
the  oxidation  reaction  on  Pt/Ru-C. 


Values  of  ii.eco  were  determined  at  the  point  of  inflection  of  the  anode  overpotential 
curves  (i.e.  1  =  ii.eco  at  d"T|a/di  =  0)  in  Region  II  for  all  the  conditions  studied  except  for  500  ppm 
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CO/H2  at  90°C  and  shown  in  Table  2.  The  steady  state  fonn  of  Equations  5  and  6  along  with 
Equation  8  above  for  itotai =  E  where  solved  simultaneously  for  0co,  Oh,  and  r\a  with  these  values 
of  ii,eco  for  various  values  of  bfc  and  kfh  such  that  a  minimum  OF  was  detennined  in  Region  II. 
The  values  determined  for  ka  will  be  discussed  in  detail  in  a  following  section.  With  these 
values  of  bfC,  0co,eq  was  determined  using  Equation  1 1  and  plotted  versus  bfC  shown  in  Figure  #. 
The  Temkin  parameters,  8(AGco),  AGcol  ,  from  Equations  15  and  16,  can  now  be 

I0CO  -0 

detennined  from  a  best  fit  to  the  bfC,  0co,eq  plot  shown  in  Figure  #. 

Estimating  ktc,  kec,  and  J3C0:  Region  III  (CO  oxidation) 

In  Region  III  as  shown  on  Figure  2,  the  slope  of  the  polarization  curve  has  decreased 
substantially.  This  indicates  that  there  is  now  available  surface  for  the  TE  oxidation  reaction  to 
occur  normally.  Also  in  this  region  4  contributes  most  to  the  total  rate  ’  (4  +  ico )  because  the 
flux  of  CO  (i.e.  ic)  is  small  due  to  the  relative  concentration  of  CO  in  the  feed.  Note  that  Oco 
approaches  zero.  Even  though  the  CO  oxidation  reaction  may  not  contribute  significantly  to  the 
total  rate,  Oco  decreases  significantly  from  this  reaction. 

The  smaller  the  value  of  kfC  the  more  TE  adsorption  controls  Region  IV  i.e.  the  value  of 
kfC  must  be  chosen  such  that  ico  approaches  zero  in  that  region  (see  Region  IV  below).  Therefore 
a  criteria  of  100  pA/cm  as  the  maximum  ico  was  assumed  at  the  highest  CO  concentration  (1% 
CO  in  TE)  for  which  kfC  was  determined.  The  steady  state  fonns  of  Equations  5  and  6  and  7  and 
8  along  with  itotai  =  ico+ih,  were  solved  simultaneously  at  a  particular  current  density  for  the  anode 
overpotential  (rja),  Oco  and  Oh  assuming  a  value  of  0.75  for  (3co  (aco  =  0.5)  for  various  values  of 
kfC  and  kec  until  a  minimum  value  of  OF  was  obtained.  However  these  values  of  kfC  and  kec 
although  optimized  may  generate  ico  which  depending  on  its  size  did  not  allow  for  the  existence 
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of  a  Region  IV  as  discussed  above.  Note  that  the  values  of  kec  and  (3C0  used  in  determining  kfC  are 
not  optimized  form  the  data  however  they  were  sufficient  for  that  purposes  of  estimating  kfC. 

Now  with  a  value  of  kfc  known  the  same  set  of  equations  from  above  were  again  used  for  various 
values  of  kec  and  j8co  until  a  minimum  OF  was  obtained  for  all  conditions  except  the  ones  at  500 
and  3,000  ppm  CO/FF  at  101  kPa  and  500  ppm  CO/H2  at  90°C  and  202  kPa  due  to  the  lack  of 
points  in  Region  III. 


Estimating  kjp.  Region  IV  (FF  adsorption  limited) 

A  sharp  increase  in  the  overpotential  similar  to  that  observed  in  Region  II  can  be  seen 
from  Figure  2  in  this  region.  However  the  CO  present  on  the  surface  has  been  fully  oxidized 
( Oco  =  0)  at  this  point  meaning  that  there  should  be  enough  available  surface  to  allow  for  the  H2 
oxidation  to  proceed  at  the  rate  of  hydrogen  adsorption.  Note  that  the  intrinsic  rate  (i.e.  per  site) 
is  dependent  of  ct  and  that  the  total  rate  of  hydrogen  adsorption  me  be  limited  by  the  residual 
CO.  unrestricted.  That  is  some  of  the  CO  may  permanently  poison  some  of  the  sites.  Now  since 
the  H2  reaction  becomes  adsorption  limited,  the  H2  surface  coverage  approaches  zero  therefore 
fixing  and  itotai  to  i/iP  reducing  Equation  6  to: 


—  =  kfhPH 


at  ico  =  0 


[27] 


We  determine  ii,p  from  Equation  31  taking  the  values  of  kfh  detennined  from  Region  II.  If  these 
values  are  lower  than  the  last  current  density  value  in  the  particual  data  being  analyzed  then  bfc 
will  be  re-estimated.  All  the  values  for  ii,eco  and  i/iP  determined  using  the  above  analysis  and  used 
in  the  parameter  estimation  are  shown  in  Table  #  below. 


Anode  Overpotential  and  CO  Isotherm  Comparisons 


13 


The  anode  overpotential  rja  is  determined  by  solving  the  steady  state  from  of  Equations  5 


and  6  including  the  current  expressions  (Equations  7  and  8)  with  the  parameters  estimated 
according  to  the  analysis  shown  above  with  values  listed  in  Tables  1  and  2.  Figure  3  a  and  b 
show  the  reliability  of  the  model  and  verifies  the  assumption  pertaining  to  the  H2  oxidation 
reaction  mechanisim.  However  since  they  were  detennined  using  the  ilp  values  from  the  data, 
we  are  unable  to  see  the  effect  of  the  isothenn. 

In  order  to  identify  the  best  isotherm  model  for  the  given  points  on  Figure  3  we  start  by 
comparing  two  well  known  Langmuir  and  Temkin  isotherm  models. 

If  this  in  fact  was  a  Temkin  isothenn,  the  value  of  AE  will  not  change  as  the  coverage 
changes.  Therefore  to  check  for  this,  the  previously  determined  values  of  CO  coverage  in  Table 
I  were  used  solving  Equation  28  for  bfc,  and  solving  this  and  Equation  30  for  AE.  These  values 
of  AE  shown  in  Table  II  are  not  constant  for  the  data  range  studied  which  and  change  quite 
considerably  especially  for  the  500  ppm  CO/H2  at  90°C  data  set.  The  same  is  true  in  the  case  of 
the  Langumir  model.  Which  indicate  a  different  isotherm  may  be  driving  the  process. 


Figures  5a  and  b  and  6a  and  b  show  a  comparison  in  the  anode  overpotential  predictions 
using  the  bfc  values  from  Table  along  with  the  ii>p  values  from  Table  I  with  those  using  the 
Langmuir  model  (Equation  28)  for  70  and  90°C  respectively.  The  predictions  are  not  accurate  in 
Region  I  where  the  anode  overpotential  is  the  lowest  and  is  expected  since  the  error  in  those 
points  as  a  result  of  the  method  used  in  determining  them  in  that  region  is  usually  pronounced. 
However  they  fit  reasonably  in  the  Regions  II  -  IV,  which  indicates  that  model  and  the 
parameters  used  that  were  estimated  in  this  study  are  adequate  for  capturing  the  phenomena. 
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Observe  that  the  curves  using  the  Langmuir  model  generally  fit  the  curves  close  to  those  using 
Equation  27  at  the  higher  partial  pressures  in  Region  II,  which  can  be  expected  upon  observation 
of  Figure  4.  However  there  is  a  large  deviation  at  the  lower  partial  pressures  in  that  region 
especially  for  500  ppm  CO/H2  at  70°C  and  is  expected  since  the  Langmir  model  was  unable  to 
predict  that  point  on  the  isothenn  graph  (Figure  4).  Also  note  the  severe  deviation  in  the 
predicted  anode  polarization  in  the  case  of  500  ppm  CO/H2  at  90°C  from  that  of  the  actual 
polarization.  This  may  be  because  of  internal  air-bleeding  which  is  a  result  of  additional 
however  non-faradic  CO  oxidation  by  the  O2  crossing  over  from  the  cathode  side.  All  the 
predictions  did  not  show  significant  differences  for  Regions  III  and  IV  and  is  expected  since  the 
pennanent  poison  phenomena  is  well  defined. 

The  anode  overpotential  curves  for  the  1  atm  data  sets  3,000  and  500  ppm  CO/H2 
collected  in  Ref.  1  along  with  the  model  predictions  using  the  parameters  in  Table  2  are  shown  in 
Figure  #.  The  values  of  ilco  were  detennined  directly  from  the  data  using  a  previously  mention 
analysis  and  ilp  was  determined  based  on  the  isothenn  in  Figure  #.  Upon  observing  Figure  #  of 
Ref.  1  we  can  see 

Poisoning  and  Recovery  Comparisons 

Figures  8-9  show  the  predictions  of  the  model  using  the  parameters  estimated  in  this 
work  shown  in  Table  2.  These  figure  used  the  Ea(js  associated  with  the  actual  CO  coverage  used 
in  constructing  Figure  3  and  the  permanent  poisoning  parameters  according  to  Figure  6  since 
they  gave  the  best  predictions  for  the  steady  state  anodic  polarization  curves.  All  the  predictions 
seem  to  follow  the  general  trend  of  each  curve  although  however  they  seem  to  under  predict  and 
over  estimate  the  cell  voltage  in  some  areas.  This  is  mainly  attributed  to  the  fact  that  there  are 
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some  differences  between  the  cell  voltage  from  the  steady  state  anode  polarization  curves  and 
that  determined  during  poisoning  and  recovery  for  most  of  the  data  set  and  not  because  of  the 
parameters.  It  is  important  to  note  that  final  cell  voltage  values  were  different  for  each  CO 
partial  pressure  investigated  and  did  not  vary  in  a  predictable  way  for  the  various  runs  at  the  end 
of  recovery.  The  predictions  were  unable  to  accurately  predict  these  final  cell  voltage  values. 
This  may  be  because  a  fixed  value  was  used  for  the  total  density  of  sites,  which  may  be  different 
for  each  run. 

The  volume  of  the  humidity  bottle  used  for  the  experiments  in  Ref.  8  was  measured  to  be 
750  ml.  The  flow  rate  into  the  bottle  was  calculated,  based  on  the  analysis  presented  by  W-k  Lee 
(21),  to  be  1.68  cm3/sec  at  70°C. 


Pco  =  Pco  (=0  exP 


V  TB  J 


[##] 


Therefore  the  maximum  residence  time  is  around  446  seconds  considering  that  the  full 
volume  of  the  bottle  is  utilized.  As  can  be  seen  from  the  graph  60%  of  this  value  (200  seconds) 
gave  the  best  fit  to  the  curve  validating  the  assumption  that  the  partial  pressure  of  CO  in  the 
humidity  bottle  can  be  acts  as  a  time  dependent  CSTR  and  that  the  humidity  bottle  is  where  the 
gas  spends  the  most  time. 


Conclusions 

The  Heyrovsky-Volmer  mechanism  was  shown  to  be  applicable  for  H2  oxidation  on  a 
Pt/Ru-C  catalyst.  A  Langmuir  model  was  adequate  for  capturing  the  CO  coverage  for  the  range 
of  partial  pressure  and  temperatures  studied  particularly  at  higher  partial  pressures,  which  is 
typical  of  this  type  of  isothenn.  Region  IV  played  an  important  role  in  the  parameter  estimation 
process  as  it  revealed  the  effect  of  permanent  poisoning  of  the  overall  adsorption  rate  constant 
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for  H2  adsorption  necessary  for  predicting  the  limiting  currents  observed  during  CO  oxidation. 
The  sensitivity  of  the  model  to  the  isothenn  validates  that  the  adsorption  and  oxidation  process 
of  H2  and  CO  are  enough  to  capture  the  effect  of  CO  poisoning  in  a  PEMFC.  The  poisoning  and 
recovery  process  could  be  modeled  according  to  the  development  in  the  text  with  these 
parameters. 
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Cell  voltage  (V) 


2 

Figure  1.  Transient  performance  with  neat  hydrogen  and  3000  ppm  CO  at  600  mA/cm  with 
CARBEL™  CL  GDM.  (TceU  =  70°C  and  P(A/C)  =15/15  psig) 
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Anode  Overpotential  [Vcell-Vcellco/H2]  (V) 


Figure  2.  Regions  identification  graph  (Data  from  Figure  3  of  Ref.  1  for  3000  ppm  CO  in  H2  at 

70°C  and  202  kPa). 
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Figure  10.  bfc  versus  0co  for  500,  3,000  and  10,000ppm  of  H2/CO  feed  at  70  and  90°C 
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Anode  Overpotential  [Vcell-Vcellco/H2]  (V) 


Figure  3a.  Overpotential  model  comparison  for  (A)  500,  ( I  )  3,000  and  (  •  )  10,000  ppm 
H2/CO  at  a  cell  temperature  of  70°C  (Lines  are  computed  from  Equations  5-8  using  bfc  from 

Table  I  and  parameters  from  Table  III). 
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Anode  Overpotential  [Vcell-Vcellco/H2]  (V) 


Figure  3b.  Overpotential  model  comparison  for  (A)  500,  ( I  )  3,000  and  (  •  )  10,000  ppm 
H2/CO  at  a  cell  temperature  of  90°C  (Lines  are  computed  from  Equations  5-8  using  bfc  from 

Table  I  and  parameters  from  Table  III). 


23 


Figure  4.  CO  Isotherm  comparison  graph  for  CO  for  500,  3,000  and  10,000ppm  of  H2/CO  feed 
at  (  •  )  70  and  (A)  90°C  (Temkin  model  (dashed  line)  and  Langmuir  model  (line)). 
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Anode  Overpotential  [Vcell-Vcellco/H2]  (V) 


Figure  5a.  Overpotential  model  comparison  for  (A)  500,  ( I  )  3,000  and  (  •  )  10,000  ppm 
H2/CO  at  a  cell  temperature  of  70°C  (Lines  are  computed  from  Equations  5-8  using  Lamgmuir 
Model  i.e.  bfc  from  Equation  12  and  parameters  from  Table  III). 
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Anode  Overpotential  [Vcell-Vcellco/H2]  (V) 


Figure  5b.  Overpotential  model  comparison  for  (A)  500,  ( I  )  3,000  and  (  •  )  10,000  ppm 
H2/CO  at  a  cell  temperature  of  90°C  (Lines  are  computed  from  Equations  5-8  using  Lamgmuir 
Model  i.e.  bfc  from  Equation  12  and  parameters  from  Table  III). 
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Anode  Overpotential  [Vcell-Vcellco/H2]  (V) 


Figure  5a.  Overpotential  model  comparison  for  (A)  500,  ( I  )  3,000  and  (  •  )  10,000  ppm 
FL/CO  at  a  cell  temperature  of  70°C  (Lines  are  computed  from  Equations  5-8  using  Temkin 
Model  i.e.  bfc  from  Equation  14  and  parameters  from  Table  III). 
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Anode  Overpotential  [Vcell-Vcellco/H2]  (V) 


Figure  5b.  Overpotential  model  comparison  for  (A)  500,  ( I  )  3,000  and  (  •  )  10,000  ppm 
FL/CO  at  a  cell  temperature  of  90°C  (Lines  are  computed  from  Equations  5-8  using  Temkin 
Model  i.e.  bfc  from  Equation  14  and  parameters  from  Table  III). 
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Anode  Overpotential  [Vcell-Vcellco/H2]  (V) 


Figure  6a.  Overpotential  model  comparison  for  (A)  500  and  ( I  )  3,000  ppm  H2/CO  at  a  cell 
temperature  of  70°C  (Lines  are  computed  from  Equations  5-8  using  bfc  from  Table  I  and 

parameters  from  Table  III). 


29 


Table  1.  Summary  of  Data  Used  for  Regional  Analysis 


Temperature 

(°C) 

Pressure 

(kPa) 

CO  Concentrations 

(ppm  CO/H2) 

70 

101 

500 

70 

101 

3,000 

70 

202 

500 

70 

202 

3,000 

70 

202 

10,000 

90 

202 

500 

90 

202 

3,000 

90 

202 

10,000 

Table  2.  Summary  of  Parameters  Determined  from  Regional  Analysis 


Parameters 

Region  I: 

bfh,  keh 

Region  II: 

bfc(S(AGC0),  AGco|eco=0,A) 

Region  III: 

kfc,  keh,  Pco 

Region  IV : 

kfh 

Table  1.  Value  of  ii,p,  ii,eco  and  bfC  determined  from  data 


CO 

Concentration 

(ppm) 

Temp  (C) 

Pco 

(Pa) 

Il,p  2 
(A/cm2) 

ii.eco 

(A/cm2) 

bfc/Pco 

500 

70 

85.5 

2 

0.328 

0.196302 

3,000 

70 

513 

1.65 

0.088 

0.056836 

10,000 

70 

1710 

1.15 

0.0399 

0.035943 

30 


500 

90 

72 

2.1 

N/A 

20 

3,000 

90 

432 

1.7 

0.402 

0.309509 

10,000 

90 

1440 

1.3 

0.152 

0.132504 

Table  2.  Comparison  of  CO  Coverage  from  Data  to  Springer  et  al.  values 


CO 

Concentration 

(ppm) 

Temp  (C) 

Pco 

(Pa) 

9co 

(Nwoga) 

0co 

(Springer) 

AE 

(J/mol) 

500 

70 

85.5 

0.835 

0.994 

879.7 

3,000 

70 

513 

0.946 

0.999 

2442.3 

10,000 

70 

1710 

0.965 

0.999 

4598.0 

500 

90 

72 

0.047 

0.969 

206819.3 

3,000 

90 

432 

0.763 

0.994 

3499.8 

10,000 

90 

1440 

0.883 

0.998 

4242.6 

Table  3.  List  of  Parameters  Used  in  Model 


bfh/Pm 

4 

kec0~ 

5.16  x  10'5  ±  (5  x  10'3)  mol  cm'2  s'1 

Eads 

72.6  ±  (5  x  lO'^kJmol'1 

Ea,kec 

60.8  kJ/mol 

A 

2.95  x  10'13  Pa1 

fiii 

0.75 

kfc 

0.2  x  10'b  ±  (5  x  10'3)  mol  cm'2  s'1  atm'1 

fico 

0.6  ±  (5  x  10'3) 

i  c;eam 

Kfh 

7.34  x  10'b  ±  (5  x  10'3)  mol  cm'2  s'1  atm'1 

keh 

1.3  x  10'5  ±  (5  x  10'3)  mol  cm'2  s'1 

Figure  Captions 

Figure  1.  Transient  performance  with  neat  hydrogen  and  3,000  ppm  CO  at  600  mA/cnr  with 
CARBEL™  CL  GDM.  (Tceu  =  70°C  and  P(A/C)  =15/15  psig). 


Figure  2.  Regions  identification  graph  (Data  from  Figure  3  of  Ref.  1  for  3000  ppm  CO  in  FE  at 
70°C  and  202  kPa). 
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Figure  3.  CO  Isotherm  comparison  graph  for  CO  for  500,  3,000  and  10,000ppm  of  H2/CO  feed  at 
(  •  )  70  and  (A)  90°C  (The  lines  are  computed  from  Equations  28  and  29  with  the  parameters  in 
Table  III). 

Figure  4.  Graph  for  values  of  iiiP/ii,pcIean  as  a  function  of  CO  partial  pressure  at  70  and  90°C.  Data 
points  ( I  )  and  (  •  )  are  for  70  and  90°C  respectively. 

Figure  5a.  Overpotential  model  comparison  for  (A)  500,  (  I  )  3,000  and  (  •  )  10,000  ppm 
FE/CO  at  a  cell  temperature  of  70°C  (Lines  are  computed  from  Equations  5-8  using  bfc  from 
Table  I  and  parameters  from  Table  III). 

Figure  5b.  Overpotential  model  comparison  for  (A)  500,  ( I  )  3,000  and  (  •  )  10,000  ppm 
FE/CO  at  a  cell  temperature  of  70°C  (Lines  are  computed  from  Equations  5-8  using  Lamgmuir 
Model  i.e.  bfc  from  Equation  29  and  parameters  from  Table  III). 

Figure  6a.  Overpotential  model  comparison  for  (A)  500,  (  I  )  3,000  and  (  •  )  10,000  ppm 
H2/CO  at  a  cell  temperature  of  90°C  (Lines  are  computed  from  Equations  5-8  using  bfc  from 
Table  I  and  parameters  from  Table  III). 

Figure  6b.  Overpotential  model  comparison  for  (A)  500,  (  I  )  3,000  and  (  •  )  10,000  ppm 
H2/CO  at  a  cell  temperature  of  90°C  (Lines  are  computed  from  Equations  5-8  using  Lamgmuir 
Model  i.e.  bfc  from  Equation  29  and  parameters  from  Table  III). 
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